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It is shown how initial conditions can be appropriately defined for the integration 

of Lorentz-Dirac equations of motion. The integration is performed forward in time. 

The theory is applied to the case of the motion of an electron in an intense laser pulse, 



has put into focus the classical description of the dynamics of relativistic electrons. Under the action of 

high intensity electromagnetic field, a major ingredient of the dynamics is the electron self-interaction 

> . . . 

JTi which implies in the damping of the movement caused by the interaction of the charge with its own 

field. The derivation of the damping force has been reviewed recently pi] revealing its relativistic origin 

asociated to the asymmetry introduced by the Doppler effect in the forward and backward emission 

of radiation. The inclusion of this force in the equation of motion leads to the nonlinear covariant 

. y ■ Lorentz-Dirac(LD) equation j|] for a point charge. 

£^ Nonlinear effects have been observed in recent experiments with intense-laser relativistic-electron 

P-i- 

• . scattering at laser frequencies and field strengths where radiation reaction forces begin to become 

significant ^,Q. Relativistic nonlinear Thompson scattering has also been observed ||. These ex- 

h : 

periments justify the recent attention the classical electrodynamics theory(CED) has received. A 
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self-consistent classical theoretical treatment of the radiative reaction force would also be useful in 
simulating future electron accelerators 0. Further, the kind of study reported upon in this paper 
may be useful in the quantal treatment (71 . 

A certain number of problems, conceptual and practical, is known to be associated with LD equation. 
These difficulties may be traced to the fact that it contains a dependence on the derivative of the 
acceleration which implies in the necessity of imposing, in order to solve it, an extra condition besides 
the usual initial conditions on position and velocity of classical mechanics. It has been estabilished 
by Rohrlich that this condition is given by the asymptotic constraint that at the far future, when 
the fields vanish, the acceleration should also vanish. The solution which is obtained however, when 
extrapolated to the instant when the external force is applied violates causality. On the other hand, 
from the practical point of view, a condition put in the far future is awkward to be implemented 
specially in a scattering situation. 

A recent book || has been devoted to this question of causality violation. A clear explanation to 
its origin was given and it has been shown how it can be solved in the context of the CED. On the 
other side, attempts has been done in order to replace Rohrlich condition by an equivalent initial 
condition PjlQ]. The purpose of this letter is to apply these modern advances, to the case of the 
classical description based in LD equation, of the movement of an electron interacting with a short, 
strong laser pulse. This problem has been recently discussed in the literature Ref. B, which provided 
numerical solutions of the Lorcntz-Dirac equation. The integration was performed backwards in time 
so that the unphysical, exponentially growing homogeneous solutions of LD would damp out, resulting 
in a numerical stable solution. 

We are going to show that Lorentz-Dirac equation of motion can be integrated forward in time with 
conditions specified at t = 0. The idea is to construct the series solution of LD equation. The initial 
acceleration is then provided by replacing in the series the velocity by their initial value at the instant 
when the external force is applied. It is easy to show that the solution obtained with this procedure, 
when extrapolated to the distant future, satisfies Rohrlich condition. However, we still have to cope 
with the existence of the unphysical runaway solutions which although formally eliminated troubles 



the process of numerical integration. Then, by combining the recursive use of the series solution with 
implicit methods of numerical integration we show that the process of integration forward in time can 
be performed. 

We write Lorentz-Dirac (LD) equation of motion as 

fd 2 v 



F l ? - a »- e [l^ +W X ) C 1 ) 

where v^^a^ and F^ xt are, respectively, the four- vector components of the velocityacceleration and of 
the external force given explicitly by 

»„=7(l,/3), (2) 

1^=7(1,-/3), (3) 






and 



F-*= 7 (^.F ext ,F ext ). (4) 

In these equations, r is the dimensionless proper time dr — LUydt/j, 7 is the relativistic factor 7 = 
. 1 with (3 = - 4| and wo is the frequency of the laser pulse with which the electron is interacting. 
We follow here the same units of Ref. ||. With these definitions, it can be easily verified that v^ L a^ = 
as it should. The quantity e = wo7"o, where tq — e 2 /m e c :i = 0.626 x 10 _23 s is the Compton time scale. 
We are assuming that the external force applies at r = and the second term in the right hand side 
of the equation represents the damping force arising as the charge starts to radiate pT| . 

By contracting the Lorentz-Dirac equation with a^, it is found that (a^a 11 ) (r) = Cexp(2r/e) is a 
solution of the resulting homogeneous equation (no external force) for times greater than To- These are 
the so-called runaway solutions. They are eliminated by imposing the Dirac-Rohrlich (DR) condition 



r — >oo 



lim a^ = (5) 

when 



lim F^ (r) = 0. (6) 

In contrast to the usual initial value problem encountered in classical physics, where all quantities, 
position, velocity, are fixed at t = 0, the above asymptotic condition, Eq. (^|) recast the problem 
into a boundary value one where the x,v are known at r = and a is forced to be zero at r = 
oo, in accordance with the condition on the force, Eq. (|6|). A particular situation is that of Ref. 
|g|, in which the problem was solved in the rest frame of the particle at some future time with 
three final homogeneous conditions. In both situations however, the numerical integration has to be 
performed backwards in time. Although this procedure eliminates the unphysical runway solutions, it 
is uncomfortable, from the point of view of applications, to have a condition given at some final time. 
It is therefore desirable to have an equivalent condition on the acceleration defined at the initial time. 

For an electron moving in an electromagnetic field the external force is given by 

FJ l xt = -{d v A li -d^A v )v v (7) 

where the quantity A^ is the vector potential given in units of moc/e. For a linearly polarized laser 
pulse, 

A„ = ($/c,A), A = xA x {<t>), $ = (8) 

which is a function of the invariant phase of the traveling wave, 

4> = Wx^t) =x -z , (9) 

where /c M = (1, 0, 0, 1) is the dimensionless laser wave number. 

Following Ref. g] we use <j> as the independent variable to recast the Dirac Lorentz equation as 

d 2 v x ( (dv\ 2 fd-f\ 2 \ 1 fdv x (\ du\ A ,\ 
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where the laser pulse electric field 

dA x {<j>) 



AG{<j>)sin(j) :, (12) 



has been introduced, with A being the maximum amplitude of the pulse, G((f>) = e (^/ A< W is a unit 
Gaussian envelope of width A(f>. Note that 
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and 
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and 

du dj dv z 
d(j) d<j) d(j) 

Substituting these relations into Eqs. (p] ) and (O) we find 

d 2 v x 1 (\dv x v x (dv x \ 2 ( v z \ dv x dv z A . \ 

d 2 v z 1/1 dv z v x dv x dv z ( v z \ I dv z \ v x A 
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7 = VI + ^ + « 2 2 (19) 

M = 7 - u z . (20) 

In Ref. H, by specifying final homogeneous conditions on the acceleration and the velocity and 
then integrating backward in time, the solution to these equations were obtained at all times. We now 
want to show that this problem can also be solved by specifying initial conditions on the motion and 
integrating forward in time. 

For the "initial" velocity of our method, we use the final velocity of the backward integration method 
of @] . As for the initial acceleration we employ the first terms of the series generated by expanding 
the equations of motion, in terms of the small quantity s = luoTq. To obtain this series we write the 
two components of the equations of motion (^6|) and (|l7]) as 
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From the above equations we derive the zeroth order for the derivatives of the components of the 

acceleration 

/,#„, \ d 
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Substituting these relations back into (Elh and ( p2| ) , we obtain the first order approximation for the 

components of the acceleration ( £& J , I ^- J . On the other hand, the zeroth order of the second 

derivative of the components of the acceleration are given by 
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Taking now the derivatives of (pTh and d22) we find the rather lengthy relations 
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where Q' , 7' and u' are the derivatives of these quantities which can be easily obtained. With these 
expressions inserted in equations (El]) and (G3) , we obtain the next order approximation for the 
acceleration. This process can be repeated to generate higher order corrections as desired. 

The forward integration using this above procedure for the initial acceleration was performed using 
the subroutine Stiffs taken from the Numerical Recipes. The efficiency of the procedure is highly 



improved when this perturbative expression of the acceleration is recurrently used at each step of the 
integration process. The results of the numerical calculations are shown in Figs. 1 and 2, respectively, 
for the transversal and longitudinal components of the momentum. The parameters are A = 1, 
e = 0.05 and Ac/) — 10. These results are identical to those obtained in Ref. j| using the backward 
integration method. 

In conclusion, it has been shown that LD equation can be integrated forward in time with an 
appropriate acceleration initial condition. This acceleration is provided by the equation of motion 
itself, treating the radiation damping term as a perturbation. This peculiar behaviour of LD equation 
stems from the fact that although the self-interaction term turns it into a third-order differential 
equation, only its particular solution is physically meaningful: all homogeneous solutions have to 
be excluded. This constraint makes the solution unique. We have discussed the forward numerical 
integration procedure in the important case of the nonlinear Compton scattering. The essential feature 
of the calculations is the crucial role played by the recurrent use of the perturbative series expansion of 
the acceleration in order to achieve stable and accurate results. We may conclude that the final Dirac- 
Rohrlich constraint can indeed be replaced by an equivalent initial condition as long as the radiation 
damping term can be treated as a local perturbation in the iterative solution of the equation of motion. 
This does not mean that nonlinear effects have been neglected. Indeed, by perturbatively iterating the 
equation of motion, a local expression for the acceleration is generated, in which the radiation affects 
the motion as much as it is affected by it. The results presented here clearly show that nonlinear 
LD equations can still be integrated forward in time if the appropriate, albeit unconventional, initial 
conditions were utilized. 
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Figure Captions: 

Fig. 1 The transversal component of the momentum vs. proper time, obtained with the forward 
integration method. See text for details. 

Fig. 2 The same for the longitudinal component of the momentum. 
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